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We calculate the reduced single-particle density matrix (RSPDM), momentum distributions, nat- 
ural orbitals and their occupancies, for a strongly interacting one-dimensional Bose-Fermi mixture 
in a double-well potential with a large central barrier. For mesoscopic systems, we find that the 
ground state properties qualitatively differ for mixtures with even number of particles (both odd-odd 
and even-even mixtures) in comparison to mixtures with odd particle numbers (odd-even and even- 
odd mixtures). For even mixtures the momentum distribution is smooth, whereas the momentum 
distribution of odd mixtures possesses distinct modulations; the differences are observed also in the 
off-diagonal correlations of the RSPDM, and in the occupancies of natural orbitals. The calcula- 
tion is based on a derived formula which enables efficient calculation of the RSPDM for mesoscopic 
mixtures in various potentials. 

PACS numbers: 05.30.-d, 03.75.Ss, 67.85.Pq 



I. INTRODUCTION 

The experiments with interacting ultracold atomic 
gases loaded in various external potentials offer great op- 
portunities for probing versatile many-body states (e.g. 
see [l| for a review), from weakly interacting gases up 
to strongly correlated states. Interactions between the 
atoms can be tuned in some cases (e.g., by employing 
Fcshbach resonances), while external potentials can as- 
sume various shapes including optical lattices, elongated 
and transversely tight traps and many others. Systems of 
interacting atoms in double-well potentials exhibit par- 
ticularly interesting phenomena which were studied over 
the years, for example, a bosonic Josephson junction 
1, 0, IE] , squeezing and entanglement of matter waves 
matter wave interference [H, Q , and recently exact 
many-body quantum dynamics in a one dimensional (ID) 
quantum well [Io| . In this work we focus on a strongly in- 
teracting Bose-Fermi mixture of particles in a ID double 
well potential. 

The achievement of quantum degeneracy in Bose- 
Fermi mixtures [TO, QJ, Gj, Q EE HE G3, H| has stim- 
ulated many studies of these systems. In ID geome- 
try, several theoretical approaches explored such mix- 
tures. For example, the mean-field a ppr oximation has 
been utilized to study phase separation [191 ]. However, for 
strongly correlated systems, which are more likely to oc- 
cur in one- than in three-dimensional systems, the mean 
field approach is not appropriate. These systems can be 
studied by using exactly solvable models and/or sophis- 
ticated numerical calculations. Luttinger liquid theory 
has been used to study pairing instabilities and phase 
diagrams [IE HH ■ Numerical calculations were used to 
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obtain phase diagram for mixture with unequal masses 
[22| . One-dimensional Bose-Fermi mixture with a finite 
coupling strength and without an external trapping po- 
tential were studied in Refs. [23|, [24], Hfl ■ Recently, an 
exactly solvable model describing ID Bose-Fermi mix- 
ture with strong interactions has been studied in Ref. 
[26[ | . The ground state wave functions for arbitrary ex- 
ternal potentials were constructed; in the model, strong 
("impenetrable core") interactions are present between 
bosons, bosons and fermions, whereas fermions are mu- 
tually noninteracting and they are spin-polarized. The 
correlation functions including the one-body density ma- 
trix were addressed for the ring geometry and the har- 
monic confinement [26j . The ground state properties and 
expansion dynamics were further explored in Ref. [27j . 

The solution of the model presented in Ref. [26| follows 
the Fermi-Bose mapping idea to calculate exact wave 
functions in the so called Tonks-Girardeau (TG) model 
of " impenetrable-core" bosons [28|] . This model has been 
experimentally realized several years ago [2^ . |3G [. with 
atoms in tight transversely confined atomic waveguides 
[3lj . at low temperatures, and with strong effective in- 
teractions [3J [H, HH . Besides the wave functions [28[ , 
the correlation functions such as the reduced single par- 
ticle density matrix (RSPDM) and related quantities in- 
cluding distributions of momenta, natural orbitals and 
their occupancies, have been studied for the TG system 
over the years [M \M, H, IE H [IlEpJp, [Z| \M, El 
for the ground states on the circle 34T"3a|. in harmonic 
confinement [3p . l37l . [H, [39|, for excited " dark-soliton" 
cigenstates [43| |. in a split-trap potential O . and also 
for time-dependent states (e.g., see [iol. |41|. l42j|). 

Here we study ground state properties of a strongly 
interacting Bose-Fermi mixture in a double-well poten- 
tial; we use the model from Ref. [26[. As a first step, we 
derive a formula which enables efficient calculation of the 
RSPDM for mesoscopic mixtures in various potentials; it 
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reduces to a calculation of RSPDM for TG bosons [42j 
in an incoherent mixed state. The formula is employed 
to calculate the RSPDM, momentum distributions, nat- 
ural orbitals and their occupancies, for ground state of 
the mixture in a double- well. For mesoscopic systems, we 
find that the ground state properties qualitatively differ 
for mixtures with even number of particles (both odd- 
odd and even-even mixtures) in comparison to mixtures 
with odd particle numbers (odd-even and even-odd mix- 
tures). For even mixtures the momentum distribution is 
smooth, whereas the momentum distribution of odd mix- 
tures possesses modulations; the differences are observed 
also in the off-diagonal correlations of the RSPDM, and 
in the occupancies of natural orbitals. 



II. THE MODEL 

We study a mixture of Ng bosons and Np spin- 
polarized fcrmions in one-dimensional geometry, in a 
ground state of an external potential V(x). Bosons and 
fcrmions experience the same external potential, and 
their masses are assumed to be approximately equal. 
This condition can be satisfied for combination of bosonic 
and fermionic isotopes of the same element, such as 
39(41) j( _4o ^ or Yb with several stable isotopes, and 
86(84) # 6 _87(85) Rh ( for detailed discussion, see [H and 
references therein). Bosons interact via a very strong re- 
pulsive contact potential, that is, their interaction is in 
the Tonks-Girardcau regime. Fcrmions are mutually non- 
interacting. Interaction between bosons and fermions is 
also a very strong repulsive contact potential. The total 
number of particles is denoted with N = Np + Np. The 
ground state for this system for finite but very strong 
interactions can be approximately written as [261 ] : 

4>(XF1 , . . . , XFN F , XB1, ■ ■ ■ , XBN B ) = 

N B N F 

sgn(x B j - x B i) \ \ sgn(x B j - x Fl ) 

l<i<j<N B j = l z=l 

4>s(xfi, ■ ■ ■ ,X Nf ,X B 1, ■ ■ ■ ,XN B )- (1) 



Here. 



ip s (xi,...,x N ) = — L= ^ {-) P (j)p jl (x 1 )...<f)p jN {xN), 
ViV - PeSjv 

(2) 

denotes a Slater determinant wave function constructed 
from the single particle wave functions <fij( x ), which are 
the N lowest single-particle eigenstates of the potential 
V(x): 



d 2 4>j(x) 
dx 2 



+ V(x)4> j (x)=E j 4> j (x), 



(3) 



j = 1, 2, . . . , N = Np + N B - In Eq. (J2J, P denotes a 
permutation from the group Sjy. For the clarity of the 
exposition, it is convenient to define the following quan- 
tities. Let us consider a subset of single particle states, 



chosen from the set {4>j{x)\j = 1, ■ • • , -/V} by crossing out 
k single-particle states; let J = {ji, J2, ■ ■ ■ , jk} denote the 
indices of the crossed states, and L = {l±, I2, ■ . ■ , lN-k} 
the indices of the remaining states (obviously J n L = 0, 
JUL = {1, . . . , N}). We define the Slater determinant 
state 



rs 



(Xl, . .., XN-k) 



1 



y/(N-k)l 

X! (-) P( t>Ph( x i)--- ( t ) PiN- k ( x N-k), (4) 



Pes,, 



where P is a permutation of indices (1%, 1%, . . . , lif-k)- 
Thus, the indices upon ^w'i>i2, •••,«) denote the crossed 
out states, rather than the ones used in the Slater deter- 
minant. 

Let ipf^ 2 ''"' jk ^ = Aip^ 1 ^ 2 '"' ' jk ' denote a symmetric 
Tonks-Girardeau state obtained by acting with a unit 
antisymmetric function A = Ili<i<j<Ar-fc s S ri ( x j — x i) 

/,(jlij2,--,Jfc) 



upon ip { s Jlj2 '-' jk) . The RSPDM of the state tp, 



will be denoted by pf,Q 



TG 



pfy h '-' jk \x, y) = (N -k) J dx 2 ... <I.r x , 

x^ j2 -' jk \x,x 2 ...,x N . k )* 
x tp^ 32 ''"' : ' k \y,X2 . . .,x N - k ). 



(5) 



The quantities ips, fpTG, and Ptg etc. will refer to states 
and the correlation functions obtained from the full set 
of single particle states {4>j(x)\j = 1, . . . , N}. 



III. FORMULA FOR THE RSPDM 

We are interested in properties of the ground state of 
the Bose- Fermi mixture described by the state (JXJ) . To 
explore the one-particle observables of this Bose-Fermi 
mixture, we need to construct the RSPDM of the bosonic 
and the fermionic subsystems, respectively. The one- 
body density matrix for the bosonic part of the mixture 
is defined as 

m B ,N F (x, y) = N B J dxpi . . . dxpN F dx B 2 ■ ■ ■ dx B N B 

X 7p*(x F i, . . . ,X FNf ,X,XB2, ■ ■ ■ ,XBN B ) 

X 1p{Xpi, ■ ■ ■ 1 XFN F1 y,X B 2, ■ ■ ■ ,XBN B )- ( 6 ) 

It is straightforward to verify that T]N B ,N F {x,y) can 
be calculated from the RSPDM describing TV Tonks- 
Girardeau bosons, r)N Bt N F (x,y) — Nb/N pTG{x,y)- This 
correlation function can be efficiently calculated by using 
the procedure derived in Ref. (42| . 

The calculation of the density matrix describing the 
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fcrmionic component of the mixture, 
, NB , Np (*,v)=N F J dXF2 ... dXFNpdxm ...dx 



BN E 



is much more involved. The derivation can be reduced to 
calculations of the RSPDM for an incoherent mixed TG 
state as follows: 



X lp*(x,X F2 , ■ ■ ■ ,%FN P ,ZB1, ■ ■ ■ ,X B N B ) 

X lp(y,X F2 , ■ ■ ■ ,X FNf ,X B 1, ■ ■ ■ ,XBN B ): ( 7 ) 



N F f 

Hn b ,N F {x,y) = jrj- / dx F2 ■ ■ ■ dx FNp dx B i ■ ■ ■ dx B N E 



^ sgn(a; - x B i)sgn(y - x Bl ) 



4>t{x) 

M X F2) 

M X FN F ) 
(t>l{xBl) 



4>* N {x) 

4>* N { X F2) 

4>* n {x FNf ) 
4>*n( x bi) 



My) 

$1 {x F2 ) 



M x fn f ) 
M x bi) 



4>n{v) 

<I>n{x F2 ) 



4>n(x FNf ) 
4>n{x B i) 



M x bn b ) ■■■ 4>*n( x bn b ) M x bn b ) ■■■ 4>n{xbn b ) 

I 



(8) 



The determinants above can be expanded along their 2nd 
row according to the Laplace formula, after which the 



J 



N F f 

Hn b ,n f (x, y) = —— j dx F3 . . . dx F N F dx B i ■ ■ ■ dx B N B 

N B N 

x []sgn(x - x Bi )sgn(y - x Bi ) £ {- f +1 {- f +l 

i=l h l = l 



integral over x F2 is trivially performed: 



4>*AxF2)(j>i(x F 2)dx F 2 



M x ) ■■ 

<f>t(xF3) ■■ 


■ 0*3 -l(x) 

■ <f>*j-i{xFz) 


^*+l(- T ) 

<j>* j+ i(x F3 ) ■■ 


<P*n( x ) 

■ </>*n{ x F3) 


^(xfnf) ■■ 

4>*l(xBl) ■■ 


■ 4>*j-x{xFN F ) 

■ 4>*M x Bl) 


4>*j + l{xFN F ) ■■ 
<f>* j+ l(x B l) ■■ 


■ M( X FN F 

■ M( x Bl) 


My) ■■ 

M x Fz) ■■ 


■ &j-l( X BNb) 

4>i-i(y) 

■ (t>lM x F3) 


(t>*j + l(xBN B ) - 

<f>i+i(y) 

Ml(xF3) ■■ 


■ 4 > n( x BN b 

M(y) 

4>N(XF3) 


k(xFN F ) ■■ 

M x bi) ■■ 


■ Ml(x F N F ) 

■ Ml( x Bl) 


Mi(xfn f ) ■■ 

Mi(xbi) ■■ 


4>n(xfn f ) 
4>n( x bi) 


h{XBN B ) ■■ 


■ Mi( x BNb) 


Mi(xbNb) ■■ 


4>n(x B N b ) 



N 



N(Nf 



-[)Y.t 1 N ) B,N F M x >y)> ( 9 ) 

' 3 = 1 



where we have used J (j)*(x F2 )(f>i(x F2 )dx F2 = Sji. The (f>j has been crossed out from the Slater determinant 

index j in n%] N _ t means that the single-particle state used in the fonnula for the ground state, see Sec. M 
B '' F Thus, we have derived a recursion formula which reduces 
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the calculation of the fcrmionic RSPDM fiN B ,N F (x,y), 
to the calculation of TV fermionic correlation functions, 
^N B .N F -i{x^y) (with one fermion less in the mixture). 
By successively applying the recursive formula above, it 
is straightforward to obtain 



N F \(N B + 1)\ 
TV! 



x X! 1 (^y) 

l<3i<...<3jsr p -i<JV 



N F \N B l 
TV! 



Etfi, 
Ptg 



(jlr-iMf-l) 



(x,y). (10) 



l<h<~ <jN p -l<N 



The correlation functions 



,JN F 



0, 



/' V; ,.i " ~ ( x iV) correspond 
to a mixture with one fermion and Nb bosons; it is 
straightforward to see from the definition of (j,n b ,n f that 

O'l. — Jatf-i) (ill- J«" p -l) 



%B,1 



Ptg 



i.e., the system with just 



one extra fermion has identical RSPDM (up to a normal- 
ization constant) to the system with Nb + 1 TG bosons 
placed in the proper single particle orbitals. 

Thus, the density matrix ^n b ,n f {x, y) is equal (up to a 
proportionality constant) to a sum of density matrices of 
TG states from an ensemble. Each TG state from the en- 
semble describe Nb + 1 TG bosons; these states arc con- 
structed by choosing Nb + 1 orbitals from the full set of 
single particle states {4>j{x)\j = 1, . . . , N}. Apparently, 



there are 



N 



N E 



such states, i.e., there are 



N 



terms 



in the sum (|10p . Thus, the density matrix fJ-N B ,N F {%,y) 
is equivalent to the density matrix of Nb + 1 bosons in a 
mixed TG state; the mixed state is an incoherent super- 
position of the ground state and many excited TG states, 
each of which is constructed by some choice of Nb + 1 
orbitals as stated above. The calculation thus reduces 
to applying the algorithm of Ref. 42|. Numerical cal- 



culation becomes too time consuming if the number of 
terms in the sum (fTOj) is too large. Nevertheless, it can 
be performed efficiently for mesoscopic systems. 

From the RSPDM, one can extract observables such 
as the momentum distribution, and important quantities 
like the natural orbitals (NOs) and their occupancies. For 
example, the fermionic momentum distribution is given 
by 



1 

2^ 



dxdye 



ik(x—y) 



VN B ,N F (x,y); 



(11) 



the eigenfunctions of the fermionic RSPDM, $^(2;), are 
called the natural orbitals (NOs), 



dx/j, NBi N F (x,y)^^ 



1,2,. 



(12) 

the eigenvalues A^j are the occupancies of these orbitals. 
The bosonic momentum distribution n„(k), NOs <S> V i(x), 
and occupancies A^j, are defined by using equivalent re- 
lations for the bosonic RSPDM. 
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FIG. 1: (color online) Momentum distributions n M (fc) for 
Nf = 6 (all curves) and Nb = 0, 1, 2, 3, 4; curves are ordered 
from bottom to top and shifted by a constant for better visi- 
bility. 



IV. MIXTURE IN A SPLIT TRAP 

In this section we apply the presented formalism to 
study the RSPDM, momentum distribution, natural or- 
bitals and their occupancies for a Bose-Fermi mixture in 
a double well potential of the form 



V(x) = V ho (x)+V G (x) 



2 ~ 2 ' Voe-W'r 



V X 



(13) 



where the parameter a (Vb) denotes the width (height, 
respectively) of the Gaussian barrier which splits the har- 
monic potential. We note in passing that the presented 
result may depend on the shape of the double- well poten- 
tial (e.g., for the split-box potential) as will be discussed 
below. Here we work in dimensionless units; the con- 
nection to physical units can be made straightforwardly: 
For example, if m is the mass of bosons (which is equal 
to the mass of fermions), and x un i t the unit of space 
(which can be chosen at will), then the unit of energy is 
E un it = h 2 / (2mx 2 umt ). 

In what follows we present result of numerical sim- 
ulations for the double-well parameters a = 0.3 and 
Vb = 230; the harmonic trap frequency parameter is 
v 2 = 7.5. It should be emphasized that we have fo- 
cused our attention to the cases when the splitting po- 
tential is sufficiently high, say, when Vb is at least several 
times larger than the energy of the JVth single-particle 
state of the potential Vho(x). Thus, all our conclusions 
should be understood to hold in this limit. The momen- 
tum distribution of the bosonic component was shown 
to be equivalent to that of TV TG bosons placed within 
V(x) and this case has already been addressed in Ref. 
|44| . Thus, we first turn our attention to the behavior 
of the fermionic momentum distribution in the mixture 
n M (fc), in dependence of the number of particles. Fig- 
ure [1] shows n^(k) for TVf = 6 and Nb = 0,1,2,3,4, 
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— Nf=7,Nb=0 

Nf=7,Nb=1 

Nf=7,Nb=2 

— Nf=7,Nb=3 




Nf=7,Nb=0 

Nf=7,Nb=1 

Nf=7,Nb=2 
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FIG. 2: (color online) Momentum distributions n M (fc) for 
Nf = 7 (all curves) and Nb = 0,1,2,3: curves are ordered 
from bottom to top and shifted by a constant for better visi- 
bility. 



FIG. 4: (color online) The total momentum distributions 
n M (fc) + n v (k) for Nf = 7 (all curves) and Nb = 0,1,2,3; 
curves are ordered from bottom to top and shifted by a con- 
stant amount for better visibility. 



— Nf=6,Nb=0 
Nf=6,Nb=1 




FIG. 3: (color online) The total momentum distributions 
+ nr,(k) for Nf = 6 (all curves) and Nb = 0, 1, 2, 3, 4; 
curves are ordered from bottom to top and shifted by a con- 
stant for better visibility. 
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FIG. 5: (color online) Contour plot of the fermionic RSPDM 
for N F = 6 and N B = 1. 



while Fig. [5] shows the same quantity for Nf = 7 and 
Nb = 0,1,2,3. From these figures we observe a quali- 
tative difference in behavior of the fermionic momentum 
distribution in dependence on the parity of the total num- 
ber of particles N = Np + Nb'- If N is even, then n fl (k) 
has a smooth bell-shaped profile. In contrast, when N is 
odd, then there are nonncgligiblc modulations on top of 
the bell-shaped profile of n M (fc); in our simulations, for 
the parameters presented here, we find that the number 
of peaks in n^(k) for odd N is the same as the number of 
fermions in the mixture Nf plus two additional humps on 
the side bands of the distribution. This parity dependent 
behavior is reflected onto the behavior of the total mo- 
mentum distributions [n^(A;)+n^(fc)] which are displayed 
in Figures [3] and [¥] for the same combinations of particles 
as presented in Figs. [T] and [2] Interestingly, odd-odd 
combinations yield results very similar to the even-even 



ones, whereas even-odd combination yields results simi- 
lar to the odd-even combinations. Thus, the parity of the 
total number of particles determines the behavior of the 
fermionic component, at least for the mesoscopic num- 
bers of particles studied here. 

The behavior of the fermionic momentum distribu- 
tion results from the properties of the fermionic density 
matrix ^,^ B ^ F (x,y) which is illustrated in Fig. [S] for 
N F = 6 and N B = 1, Fig. [6] for N F = 6 and N B = 2, 
Fig. [7] for Nf = 7 and Nb = 1, and in Fig. [5] for N F = 7 
and Nb = 2. From these figures we observe that the 
most significant difference between the total even and 
odd numbers of particles occurs in the second and the 
fourth quadrant of the x-y plane: If N is even, the values 
of the fermionic RSPDM for x < < y and y < < x are 
negligible, HN B ,N F (%,y) ~ 0; In contrast to that, if TV is 
odd, there are some oscillations of RSPDM in the second 
and the fourth quadrant, in particular close to the line 
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FIG. 6: (color online) Contour plot of the fermionic RSPDM 
for N F = 6 and N B = 2. 
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FIG. 7: (color online) Contour plot of the fermionic RSPDM 
for N F = 7 and N B = 1. 
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FIG. 9: (color online) The occupancies of the fermionic 
RSPDM for the combinations Nf = 6 and Nb = 1 (squares) , 
and Nf = 6 and Nb = 2 (circles). The lines serve to guide 
the eye. 
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FIG. 10: (color online) The occupancies of the fermionic 
RSPDM for the combinations Nf = 7 and Nb = 1 (squares) , 
and Nf = 7 and Nb = 2 (circles). The lines serve to guide 
the eye. 
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FIG. 8: (color online) Contour plot of the fermionic RSPDM 
for N F = 7 and N B = 2. 



fJ-NB,N F (x, —x). These observations indicate that there is 
much greater spatial coherence between the fields at the 
two sides of the well for odd N. A similar observation 
has been made for the TG gas in a split trap [44| . 

Let us now focus on the natural orbitals. that is. their 
occupancies. The occupancies A^,, corresponding to the 
four fermionic density matrices from Figs. OH] are illus- 
trated in Figs. [9] and [10] We immediately observe that 
when the total number of particles is even, the occupan- 
cies come in pairs and they correspond to the degenerate 
natural orbitals. Namely for the total even numbers of 
particles the symmetry of the system with respect to the 
double well is naturally preserved. However, for the to- 
tal odd number of particles this is not the case, which is 
reflected in the occupancies which do not come in degen- 
erate pairs but rather decrease one by one. 

The observation on the off diagonal behavior of the 
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fcrmionic RSPDM can be underpinned analytically. To 
this end, let us assume that the trap is split by an in- 
finitely strong delta function, i.e., the the mixture is in 
the potential V(x) = V^rap + k5(x), where k — > oo; the 
potential T^rap can ^ e a harmonic oscillator trap, or it 
may have some other functional form. For the strongly- 
interacting Bosc-Fcrmi mixture in a ground state of such 
a potential we can prove the following: If the total num- 
ber of particles N is even, then [iN Bt N F (x, y) = in the 



second and the fourth quadrant of the x-y plane, that is, 
for x < < y and x > > y. 

Consecutive single-particle states in a split-trap po- 
tential are degenerate, i.e., E 2m — l — E 2m f° r 771 = 
1,2,...; moreover, degenerate eigenstates are simply re- 
lated: 2m (x) = sgn(x)0 2TO -i(x). By using this relation, 
one of the two Slater determinants ([2]) which enter the 
formula for the RSPDM ([7]) [the determinant which de- 
pends on variable x] can be written as 



0i 



( cf>i(x) sgn(x)0i(x) 
0i(x 2 ) sgn(x 2 )0i(x 2 ) 



<p N -i(x) sgn(x)0iv-i(x) ^ 
4> N -i(x 2 ) sgn(x 2 )0Ar_i(x 2 ) 



\(j)i(x N ) sgn(xAr)0i(x7v) •■• (/)n-i(xn) sgn(x N )(/) N ^i{x N ) J 



(14) 



r 



Here we have simplified the notation and labeled the co- 
ordinates as Xj, where j runs up to the total number 
of particles N, that is, we do not explicitly refer to the 
fermionic or bosonic coordinates, as it is redundant for 
the proof. Let us assume that x < < y, and that Ni 
coordinates are negative, N 2 are positive, N\ > N 2 , and 
Ni + N 2 = N. First, let us demonstrate that the de- 



terminant (|14[) is zero in this case. In order to see that, 
we assume that x, x 2 , . . . , xjy t < 0, and that the rest of 
the coordinates are positive (due to the antisymmetry of 
the Slater determinant, any other choice of positive and 
negative coordinates would yield the same result). Let 
us add the first column to the second one in (fT4"|) , then 
the third to the fourth column and so on to obtain 



J 



( 0i (a?) 
0i (x 2 ) 



0i(xj Vl ) 
0i(xiVi+i) 20i(xj Vl +i) 



0JV-l(x) 
0JV-l(x 2 ) 



0JV-l(xjVi) 

0jv-i(xat 1+ i) 20jv(xjv 1+ i) 



\ 0l(xjv) 201 (x N ) ■■■ (j) N -i(x N ) 20at_i(x7v) / 

I 



(15) 



Thus, the first Ni entries of every even column, columns all the way to the rigth; the determinant is then 
2, 4, . . . , TV, is zero. Suppose that we shift all these even proportional to 

I 



/ 0i (a?) 

01 (£ 2 ) 



03 (x) 
03(^2) 



0l(xjv 1 ) 3 (XiV 1 ) 
0l(XiV 1 +l) 03(xat 1 + i) 

\ 0i(xjv) 3 (xat) 



0jV_l(x) 
f >N~l(x 2 ) 



0jv-i(xjvJ 

i>N-l(XN!+l) 20i(xjv 1+ i) 
0jv-i(xat) 20i(xat) 

I 








20Ar-l(Xj Vl - 



20jv-i(x w ) / 



(16) 



In the upper right corner there is a block of zeros with the size iVi x N/2. Since iVi > N/2, the determinant is 
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identically zero. In fact, it is straightforward to verify by 
using the procedure outlined above that the determinant 
is exactly zero whenever N\ ^ N2, as long as x and y 
have opposite signs and N is even. On the other hand, 
if N\ = N2, then the other Slater determinant entering 
the formula for the RSPDM in Eq. [the one which 
depends on y variable], does not have equal number of 
positive and negative coordinates and therefore it is zero, 
which completes our proof. It is straightforward to see 
that the same result holds for the bosonic part of the 
mixture, and therefore for the mixture as a whole. This 
result can be interpreted as follows: If two slits were 
opened on the opposite sides of the barrier, and if the gas 
was allowed to drop from the slits, expand, and interfere, 
the interference fringes would have not been observed for 
even N and sufficiently high barrier. In fact, such an 
experiment could be used to determine the parity of the 
mixture. 

We note that this results corresponds to the observa- 
tion made in Ref. [44| on the difference in spatial coher- 
ence of the Tonks- Girardeau gas in a split-trap in depen- 
dence of the parity of the number of TG bosons. Our 
proof covers this case as well. 

Finally, it is worthy to note that the shape of the ex- 
ternal potential can influence the presented results. For 
example, for a box potential with a Gaussian as a split- 
barrier we did not find the differences between the even 
and odd numbers of particles presented above, that is, 
in all cases the momentum distribution was smooth. For 
a potential of the form oc \x\ with a splitting Gaussian 
term the differences between the odd and even numbers 
are recovered. 



V. CONCLUSION 

In conclusion, we have studied a ID Bose-Fermi mix- 
ture, where the boson-boson and boson-fermion inter- 



actions are very strong and repulsive, whereas (spin- 
polarized) fcrmions arc mutually nonintcracting; the 
atoms in each species have approximately the same mass. 
The ground state for this system for finite but very strong 
interactions in an external potential has been constructed 
in Ref. (26J. We have studied the ground state prop- 
erties of the mixture in a double-well trap, with a suf- 
ficiently high barrier. More specifically, a formula for 
the calculation of the reduced one-body density matrix 
of the fermionic (and also bosonic) component has been 
derived, which was subsequently employed to study the 
momentum distribution, natural orbitals and their occu- 
pancies. We have found that the behavior of the mo- 
mentum distribution depends on the parity of the total 
number of particles: For even mixtures the momentum 
distribution is smooth, whereas the momentum distri- 
bution of odd mixtures possesses distinct modulations. 
When the total number of particles is even, the correla- 
tions expressed by the reduced one-body density matrix 
are negligible between the two wells. 

Acknowledgments 



This work is supported by the Croatian Ministry of 
Science (Grant No. 119-0000000-1015). H.B. and D.J. 
acknowledge support from the Croatian-Israeli project 
cooperation and the Croatian National Foundation for 
Science. K.L. would like to thank Miss Antonela Ozretic 
Banovac for support. 



[1] I. Bloch, J. Dalibard, and W. Zwerger, Rev. Mod. Phys. 

80, 885 (2008). 
[2] G.J. Milburn, J. Corney, E.M. Wright, D.F. Walls, Phys. 

Rev. A 55, 4318 (1997). 
[3] A. Smerzi, S. Fantoni, S. Giovanazzi, and S.R. Shenoy, 

Phys. Rev. Lett. 79, 4950 (1997). 
[4] M. Albiez,R. Gati, J. Foiling, S. Hunsmann, M. Cris- 

tiani, and M.K. Oberthaler, Phys. Rev. Lett. 95, 010402 

(2005). 

[5] S. Levy, E. Lahoud, I. Shomroni, and J. Steinhauer, Na- 
ture (London) 449, 579 (2007). 

[6] G.-B. Jo, Y. Shin, S. Will, T.A. Pasquini, M. Saba, W. 
Ketterle, D.E. Pritchard, M. Vengalattore, and M. Pren- 
tiss Phys. Rev. Lett. 98, 030407 (2007). 

[7] J. Esteve, C. Gross, A. Weller, S. Giovanazzi, M.K. 
Oberthaler Nature (London) 455, 1216 (2008). 

[8] MR. Andrews, C.G. Townsend, H.-J. Miesner, D.S. Dur- 
fee, D.M. Kurn, and W. Ketterle, Science 275, 637 



(1997). 

[9] T. Schumm, S. Hofferberth, L.M. Andersson, S. Wilder- 
muth, S. Groth, I. Bar- Joseph, J. Schmiedmayer, and P. 
Kriiger, Nature Physics 1, 57 (2005). 

[10] K. Sa kmann, A.I. Strelt sov, O.E. Alon, and L.S. Ceder- 
baum, larXiv:0905.0902V l. 

[11] A.G. Truscot, K.E. Strecker, W.I. Alexander, G.B. Par- 
tridge, and R.G. Hulet, Science 291, 2570 (2001). 

[12] F. Schreck, L. Khaykovich, K.L. Corwin, G. Ferrari, T. 
Bourdel, J. Cubizolles, and C. Salomon, Phys. Rev. Lett. 
87, 080403 (2001). 

[13] G. Modugno, G. Roati, F. Riboli, F. Ferlaino, R. J. 
Brecha, and M. Inguscio, Phys. Rev. Lett. 297, 2240 
(2002). 

[14] Z. Hadzibabic, C.A. Stan, K. Dieckmann, S. Gupta, 
M.W. Zwierlein, A. Gorlitz and W. Ketterle, Phys. Rev. 
Lett. 88, 160401 (2002). 

[15] J. Goldwin, S. Inouye, M. L. Olsen, B. Newman, B. D. 



9 



DePaola, and D. S. Jin, Phys. Rev. A 70, 021601(R) 
(2004). 

[16] C. Silber, S. Giinther, C. Marzok, B. Deh, P.W. 
Courteille, and C. Zimmermann, Phys. Rev. Lett. 95, 
170408 (2005). 

[17] S. Ospelkaus, C. Ospelkaus, L. Humbert, K. Sengstock, 
and K. Bongs, Phys. Rev. Lett. 97, 120403 (2006); 

[18] Y. Shin, A. Schirotzek, C.H. Schunk, and W. Ketterle, 
Phys. Rev. Lett. 101, 070404 (2008). 

[19] K.K. Das, Phys. Rev. Lett. 90, 170403 (2003). 

[20] M.A. Cazalilla and A.F. Ho, Phys. Rev. Lett. 91, 150403 
(2003). 

[21] L. Mathey, D.-W. Wang, W. Hofstetter, M.D. Lukin, and 
E. Dernier, Phys. Rev. Lett. 93, 120404 (2004). 

[22] M. Rizzi and A. Imambekov, Phys. Rev. A 77, 023621 
(2008). 

[23] A. Imambekov and E. Dernier, Ann. Phys. (N.Y.) 321, 
2390 (2006). 

[24] H. Frahm and G. Palacios, Phys. Rev. A 72, 061604(R) 
(2005) 

[25] M.T. Batchelor, M. Bortz, X.W. Guan, and N. Oelkers, 

Phys. Rev. A 72, 061603(R) (2005). 
[26] M.D. Girardeau and A. Minguzzi, Phys. Rev. Lett. 99, 

230402 (2007). 

[27] B. Fang, P. Vignolo, C. Miniatura, and A. Minguzzi, 
Phys. Rev. A 79, 023623 (2009). 

[28] M. Girardeau, J. Math. Phys. 1, 516 (1960). 

[29] T. Kinoshita, T. Wenger, and D.S. Weiss, Science 305, 
1125 (2004); B. Paredes, A. Widera, V. Murg, O. Mandel, 
S. Foiling, I. Cirac, G. V. Shlyapnikov, T. W. Hansen, 
and I. Bloch, Nature (London) 429, 277 (2004). 



[30] T. Kinoshita, T. Wenger, and D.S. Weiss, Nature (Lon- 
don) 440, 900 (2006). 

[31] M. Olshanii, Phys. Rev. Lett. 81, 938 (1998). 

[32] D.S. Petrov, G.V. Shlyapnikov, and J.T.M. Walraven, 
Phys. Rev. Lett. 85 3745 (2000). 

[33] V. Dunjko, V. Lorent, and M. Olshanii, Phys. Rev. Lett. 
86 5413 (2001). 

[34] A. Lenard, J. Math. Phys. 5, 930 (1964); T.D. Schultz, 
ibid. 4, 666 (1963). 

[35] V.E. Korepin, N.M. Bogoliubov, and A.G. Izergin, Quan- 
tum Inverse Scattering Method and Correlation Func- 
tions (Cambridge, Cambridge University Press, 1993). 

[36] M.D. Girardeau, E.M. Wright, and J.M. Triscari, Phys. 
Rev. A 63, 033601 (2001); G. J. Lapeyre, M. D. Gi- 
rardeau, and E.M. Wright, ibid. 66, 023606 (2002). 

[37] T. Papenbrock, Phys. Rev. A 67, 041601(R) (2003). 

[38] P. J. Forrester, N. E. Frankel, T. M. Garoni, and N. S. 
Witte, Phys. Rev. A 67, 043607 (2003). 

[39] M. Rigol and A. Muramatsu, Phys. Rev. A 70, 031603(R) 

(2004) . 

[40] M. Rigol and A. Muramatsu, Phys. Rev. Lett. 94, 240403 

(2005) . 

[41] A. Minguzzi and D.M. Gangardt, Phys. Rev. Lett. 94, 
240404 (2005). 

[42] R. Pezer and H. Buljan, Phys. Rev. Lett. 98, 240403 
(2007). 

[43] H. Buljan, K. Lelas, R. Pezer, and M. Jablan, Phys. Rev. 

A, 76, 043609 (2007) 
[44] J. Goold and Th. Busch, Phys. Rev. A 77, 063601 (2008). 



